function output = f2(lambda)

p0 = -6;
gamma = 3;

%output = tanh((lambda^.5 - (lambda + 2*p0*gamma)^.5)/gamma) - ...
%         lambda/(9.8+gamma*(lambda^.5));
output = tanh((sqrt(lambda) - sqrt(lambda + 2*p0*gamma))/gamma) - ...
        lambda/(9.8+gamma*sqrt(lambda));
end
